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The first measurement of energy produced transverse to the beam direction at RHIC is presented. 
The mid-rapidity transverse energy density per participating nucleon rises steadily with the number 
of participants, closely paralleling the rise in charged-particle density, such that (Et) / (Nch) remains 
relatively constant as a function of centrality. The energy density calculated via Bjorken's prescrip- 
tion for the 2% most central Au+Au collisions at ^Sj^„ = 130 GeV is at least esj ~ 4.6 GeV/fm^, 
which is a factor of 1.6 larger than found at ^/s]^ = 17.2 GeV (Pb-|-Pb at CERN). 



PACS Numbers: 25.75.-q, 12.38.Mh, 13.60. Le, 13.85.Hd 



The PHENIX detector |) at RHIC, the Relativistic 
Heavy Ion Collider at Brookhaven National Laboratory, 
is designed to measure the properties of nuclear mat- 
ter at the highest temperatures and energy densities. 
For example a transition to a quark-gluon plasma has 
been predicted for energy densities on the order of a few 
GeV/fm3 [|. The spatial energy density (e) in a rela- 
tivistic collision can be estimated (following Bjorken ||]) 
by measuring the transverse energy density in rapidity, 
dET/dy, which is effectively the co-moving energy den- 
sity in a longitudinal expansion: 



dEr 
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dy To 7ri?2 



(1) 



where tq, the formation time, is usually taken as 1 fm/c, 
and ttE? is the effective area of the collision. The trans- 
verse energy {Et) is a multiparticle variable defined as: 



Et^^E, 



smt 



dET{ri)/dri = sin6'(77) dE{f])/dT], (2) 



where 6 is the polar angle, rj — — In tan 0/2 is the pseudo- 
rapidity, Ei is by convention taken as the kinetic energy 
for nucleons and the total energy for all other particles [Q , 
and the sum is taken over all particles emitted into a fixed 
solid angle for each event. Et measurements, even in 
limited apertures at mid-rapidity, provide excellent char- 
acterization of the nuclear geometry of a reaction on an 
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event-by-event basis and are sensitive to the underlying 
reaction dynamics ||^ . 

During the RHIC run in the summer of 2000, PHENIX 
accumulated close to 5 million interaction triggers for 
Au-|-Au colhsions at = 130 GeV using Zero Degree 

Calorimeters (ZDC) and Beam-Beam Counters (BBC) as 
triggering devices. The events were selected with a re- 
quirement on the collision vertex position along the beam 
axis, \z\ < 20 cm, as in the recent PHENIX publication 
on mid-rapidity multiplicity distributions j^], where fur- 
ther details are given. 

The present measurement uses a section of the elec- 
tromagnetic calorimeter (EMCal) from the PHENIX 
central-spectrometer, with front face 5.1 m from the 
beam axis. This section is part of a sampling calorime- 
ter, custom developed and built for PHENIX com- 
posed of alternating Pb and scintillator tiles (PbSc) with 
readout of individual towers, 5.54 x 5.54 cm^ in cross sec- 
tion, via wavelength shifting (WLS) fibers in a "shash- 
lik" geometry. The depth of the PbSc calorimeter is 
18 radiation lengths (Xq) which corresponds to 0.85 in- 
teraction lengths. The PbSc calorimeter has an en- 
ergy resolution of 8.2%/ y/E{GeV) ©1.9% for test beam 
electrons, with measured response proportional to in- 
cident electron energy to within ±2% over the range 
0.3 <Ee< 40.0 GeV §. 

During construction, the calibration of the calorimeter 
was set by simultaneously recording the response to laser 
excitation and to cosmic-ray muons penetrating trans- 
versely to the tower axis. The calibration was maintained 
in-situ during the run by monitoring relativistic charged 
particles from Au-|-Au collisions. The absolute energy 
scale was determined by test-beam measurements nor- 
malized to electrons with known energy. A final adjust- 
ment of the absolute energy scale was performed using 
in-situ identified electrons (p > 500 MeV/c) by shifting 
the originally measured energy/momentum (E/p) peak 
from 1.02±0.01 to 1.00. The accuracy of the absolute 
energy scale was cross-checked in-situ against both the 
minimum ionizing peak (MIP) of charged particles pen- 
etrating along the tower axis and the mass of the tt'^. 
The corrected energy distribution of EMCal clusters from 
1.0 ±0.1 GeV/c charged tracks (mostly pions) measured 
in the Drift Chamber [|l[ exhibits a clear MIP (Fig. [l^), 
as well as energy due to nuclear interactions in the ma- 
terial of the EMCal. The MIP position is in agreement 
within 2% to the value obtained in the test beam (270 
MeV). The mass of the 7r°, reconstructed from pairs of 
EMCal clusters (assumed to be photons 0) of total en- 
ergy greater than 2 GeV (Fig. |^b), is within 1.5% of the 
published value. This sets the systematic error of the 
absolute energy scale at less than 1.5%. 

The data sample for the present measurement is 
taken from the same runs used in our multiplicity mea- 
surement Q (no magnetic field), and comprises about 
140,000 events from the BBC trigger which detects [92 ± 



2{syst)]% of the nuclear interaction cross section of 7.2b 
with a background contamination of [1 ± l{syst)]% 
The transverse energy was measured using the PbSc EM- 
Cal in a fiducial aperture jryj < 0.38 in pseudorapidity 
and A(j> — 44.4° in azimuth. Et was computed for each 
event (Eq. ||) using clusters of energy greater than 20 
MeV, composed of adjacent towers with deposited energy 
of more than 3 MeV. The angle 9i is computed from the 
centroid of the cluster of energy Ei assuming a particle 
originating from the event vertex. 

The raw spectrum of measured transverse energy, 
■E^TEMC, in the fiducial aperture of the PHENIX EMCal 
for Au-|-Au collisions at y/s^ = 130 GeV is shown in 
Fig. |[ upper scale. The lower scale in Fig. || represents 
a correction of the raw -Etemc by a factor of 12.8 to cor- 
respond to the hadronic dEx / dri\ri=Q in the full azimuth. 
The 12.8 is composed of a factor of 10.6 for the fiducial 
acceptance, a factor of 1.03 for disabled calorimeter tow- 
ers and a factor, k = 1.17 ± 0.01, which is the ratio of 
the hadronic Et in the fiducial aperture to the measured 
Etemc- The k factor includes the response of the de- 
tector to charged and neutral particles emitted from the 
event vertex into the fiducial aperture, and additional 
corrections for energy in-flow from outside the fiducial 
aperture and for losses ||] . These factors were calculated 
with a GEANT §] based Monte Carlo (MC) simulation 
of the detector using HIJING as the event generator . 

For Et measurements at mid-rapidity at a collider, the 
EMCal acts as a thin but effective hadronic calorime- 
ter. Charged pions with pT < 0.35 GeV/c, kaons 
(PT < 0.64 GeV/c) and protons {pt < 0.94 GeV/c)— pr 
values which are near or above the (pt) for all 3 cases — 
stop (i.e. deposit all their kinetic energy) in the EMCal. 
For higher pT hadrons, 43% leave the MIP and 57% in- 
teract, leaving an average of ~ 65% of their energy. The 
measured -Etemc ^-"^^ ^ 0-01 of the total Et striking 
the EMCal, which is composed roughly of 40% produced 
by charged pions, 40% by photons (from 7r° and other 
decays), and 20% by all other particles (including decay 
muons). The particle composition and (pt) in HIJING 
are close to the observed values, and furthermore, the k 
factor is insensitive to reasonable variations (for instance 
varying the momenta of all particles by ±15% changes 
the overall k by less than ±2%), leading to an estimated 
systematic uncertainty in k of less than ±3% due to par- 
ticle composition and momentum. 

The main issues for the MC are the in-fiow contribution 
and losses. The losses are due to particles which originate 
within the aperture but whose decay products miss the 
EMCal (10%), or whose energy is lost due to edge effects 
(6%) or clustering (2%). The in-flow, (24 ± 1)% of the 
Et striking the EMCal, is principally of two types: (1) 
albedo from the magnet poles; (2) particles which orig- 
inate outside the aperture of the calorimeter but whose 
decay products hit the calorimeter. The in-flow compo- 
nent of k was checked by comparing the MC and the 
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measurements for events with a vertex outside the nor- 
mal range, just at and inside a pole face of the axial 
central-spectrometer magnet, 38 < z < 42 cm, for which 
the calorimeter aperture is partly shadowed. The frac- 
tion of the total energy, dE-EMc/ EemCi in bins of width 
2 towers along the z coordinate of the EMCal, zemCj is 
shown in Fig. ||a. The HIJING MC simulation agrees 
with the measured data everywhere except in the range 
-^EMC > 100 cm, which is fully shadowed by the pole, 
where the simulation shows ~20% less energy than the 
data. In Fig. the distributions of the cluster en- 
ergy. Eel, for the open aperture, zemc < —50 cm, are 
shown for both HIJING and the data and are in excel- 
lent agreement. The in- flow component of HIJING is also 
indicated as a dotted line and falls much more sharply 
than the total Ed spectrum. The residual discrepancy 
of the energy in the shadowed region, which contributes 
roughly 10% of the total signal, results in a ±(2 — 3)% 
systematic uncertainty in Et due to the uncertainty in 
the in-flow. Combining this with the uncertainty due 
to particle composition and momentum yields an overall 
factor k = [1.17 ± 0.01] ± 4% (syst), which, according to 
the MC, is independent of centrality. 

Returning to Fig. ||, the shape of the measured trans- 
verse energy spectrum shows the characteristic form of 
Et distributions in limited apertures: a peak and sharp 
drop-off at low values of Et corresponding to peripheral 
collisions with grazing impact; a broad, gently sloping 
plateau at the mid-range of impact parameters, domi- 
nated by the nuclear geometry; and then at higher values 
of Et, which correspond to the most central collisions 
where the nuclei are fully overlapped, a 'knee' leading 
to a fall-off which is very steep for large apertures and 
which becomes less steep, the smaller the aperture [ pT| . 
It should be emphasized that the correction of -Etemc 
to (IEt / drj\ri=Q by a single scale factor (predominantly 
acceptance) is valid up to the knee of the distribution, 
roughly the upper 1 percentile. Above the knee, the fall- 
off depends on the aperture and is sensitive to detector 
effects as well as statistical and dynamical fluctuations. 
Thus an actual measurement of cLEt / d'q\r,=o for A77 = 1.0 
and full azimuth would have a sharper fall-off above the 
knee. With this caveat, the uncertainty in the absolute 
energy scale (±1.5%) and the uncertainty in k of ±4% are 
combined to yield an overall uncertainty in the hadronic 
dET / di]\jj=o of ±4.5% (syst), independent of Et, where 
the statistical error is negligible. 

Mid-rapidity Et distributions are a standard method 
of defining centrality P,p]-p^ . Thus, it is important to 
determine for the present data the detailed relationship of 
transverse energy production to A^part, the number of nu- 
cleons participating in the collision (participants), which 
in earlier fixed target experiments was deduced straight- 
forwardly by measuring the energy of spectator nucleons 
and fragments in a Zero Degree Calorimeter at beam ra- 
pidity. Following a procedure used in our previous publi- 



cation on the mid-rapidity charged multiplicity {Nch) dis- 
tribution, in which a clear increase of {dNch / dri\ri=o) per 
participant with the number of participants was demon- 
strated we calculate {dET / dr]\jj^o) as a function of 
centrality in upper percentile ranges of the 7.2b Au+Au 
interaction cross section (see Table |). Figure ^ shows 
that {dET / dr]\ri=a) per participant also increases with 
Apart , closely paralleling the rise in charged particle den- 
sity (Table ||). This is better illustrated in Fig. ^ where 
the ratio {dET/dr]\^=o) / {dNch/dri\ri=a) remains constant 
at a value of ^ 0.8 GeV, independent of centrality. Com- 
parison to the measurements of WA98 |jl^ from Pb+Pb 
collisions at ^s„„ — 17.2 GeV is instructive. The 
WA98 data for mid-rapidity {dET / dr]\jnid) per partici- 
pant are shown in Fig. ^ and are essentially indepen- 
dent of A^part for A^part > 200 Q . WA98 parameterizes 
their data as dET / dr]\jnid oc A'part with a — 1.08 ± 0.06 
while the same parameterization for our data yields a — 
1.13 ±0.05. Fig. ^ also shows that (dET / drjlri^o) for cen- 
tral Au+Au coUisions at y^s^ =130 GeV is about 40% 
larger than found by WA98, yet, for both cm. energies, 
{dET / drf) / {dN ch / drj) remains constant versus centrality 
at roughly the same value, ^ 0.8 GeV (Fig. ^). 

The Bjorken energy density for Pb+Pb collisions at 
= 17.2 GeV was given by the NA49 collab- 
oration jl^. NA49 reported a value of mid-rapidity 
dET/dri\i^id = 405 GeV for the most central 2% of the 
inelastic cross section, in agreement with WA98. This 
corresponds to a value of esj — 2.9 GeV/fm'^. 
A straightforward derivation of egj from our measured 
dET/dr]\,^^o of 578+3^ GeV for the same centrality cut, 
corrected to dET / dy\y=o by a factor of 1.19 ± 0.01 from 
our HIJING MC, and taking ttR^ = 148 fm^ (i.e. R = 
1.18fmyl^/'^) gives esj = 4.6 GeV/fm^^, an increase of 
60% over the NA49 value. 

In conclusion, the mid-rapidity transverse energy den- 
sity for central Au+Au collisions, and likely the spatial 
energy density, is at least 1.6 times larger at ^s„„ = 
130 GeV (RHIC) than at = 17.2 GeV (CERN). 

The variation of the Et density per participant with 
centrality is very similar to the previously reported de- 
pendence of charged multiplicity density per participant 
at RHIC energies. These results, together with the ob- 
served constancy of {Et) I {Nch) at a value ~ 0.8 GeV, 
indicate that the additional energy density at RHIC en- 
ergies is achieved mainly by an increase in particle pro- 
duction rather than by an increase in transverse energy 
per particle. 
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FIG. 1. (a) The distribution of EMCal clusters correspond- 
ing to 1 GeV/c charged tracks (mostly pions) from Au-|-Au 
collisions, (b) The reconstructed 7r° mass from pairs of EM- 
Cal clusters with total energy > 2 GeV. 
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FIG. 2. The raw Etemc distribution measured in the 
A(/) = 44.4° azimuthal and |?7| < 0.38 polar angle fiducial 
acceptance for Au-|-Au at = 130 GeV (upper scale) 

and total hadronic dEr/driln^o (lower scale), see text. The 
solid line is the minimum bias distribution with the BBC trig- 
ger; the dashed lines correspond to the distributions for the 4 
most central bins in Table I. 
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FIG. 3. 
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(a) The fraction of -Etemc in bins of 11.08 cm 
along the EMCal Zemc direction for event vertex near a pole 
face; histogram from MC simulation, solid points from beam 
data, (b) EMCal cluster energy spectrum from HIJING MC 
(solid line), with in-fiow component (dotted line), compared 
to data (solid points). 
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FIG. 4. (a) PHENIX transverse energy density per 
participant dEx / drj\n=o / Npa,Tt for Au+Au collisions at 
ySj^„ = 130 GeV as a function of A^part, the number of 
participants, compared to data from WA98 for Pb+Pb 
collisions at ,/s„ 



= 17.2 GeV. The solid line is the N^^^^ 
best fit and the dashed lines represent the effect of the ±lcr 
Apart -dependent systematic errors for dET/drj\r,=o and Apart. 
There is an additional overall (Apart-independent) system- 
atic uncertainty of ±4.5% from dEr / dr]\r)^o and ±2.0% from 
Apart, (b) PHENIX dET/dvi\r,=o/dN^h/drt\r,^o versus Apart, 
including all systematic errors, compared to WA98. Note that 
the WA98 data in both (a) and (b) have an additional ±20% 
overall systematic error which is not shown. 



TABLE I. Average transverse energy density vs. centrality. 
The statistical errors are negligible. Errors on l^dEr / dr]\n=o) 
are the Apart-dependent systematic errors from the uncer- 
tainty of the BBC cross section such that all points move 
together. There is an additional overall (Apart-independent) 
systematic uncertainty of ±4.5%. 



Centrality {dEr / d7^\r,=o) (GeV) {dA,fa/d?7|„^o) @ (Apart) |5 | 






- 5% 


503 ±2 


622 ± 41 


347 ± 10 


5 


- 10% 


409 ±4 


498 ± 31 


293 ± 9 


10 


- 15% 


340 ±5 


413 ± 25 


248 ± 8 


15 


- 20% 


283 ±7 


344 ± 21 


211 ± 7 


20 


- 25% 


233 ±7 


287 ± 18 


177 ± 7 


25 


- 30% 


191 ±8 


235 ± 16 


146 ± 6 


30 


- 35% 


154 ±8 


188 ± 14 


122 ± 5 


35 


- 40% 


123 ±7 


147 ± 12 


99 ± 5 


40 


- 45% 


98 ±7 


115 ± 11 


82 ± 5 


45 


- 50% 


76 ±6 


89 ± 9 


68 ± 4 
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